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Abstract. For stationary, barotropic fluids in Newtonian gravity we give simple criteria on 
the equation of state and the "law of motion" which guarantee finite or infinite extent of 
the fluid region (providing a priori estimates for the corresponding stationary Newton-Euler 
system). Under more restrictive conditions, we can also exclude the presence of "hollow" 
configurations. Our main result, which does not assume axial symmetry, uses the virial 
theorem as the key ingredient and generalises a known result in the static case. In the axially 
symmetric case stronger results are obtained and examples are discussed. 



1. Introduction 

This work deals with a priori estimates of solutions of the stationary Newton-Euler system of 
equations. For definiteness the latter includes a barotropic equation of state (EOS) p g = p g (p) 
relating the matter density and the pressure, and a "law of motion" (LOM) specifying the 
velocity U(x') or the centrifugal potential (p c = C (^) as a function of position. This 
system can be used to model stars or galaxies, provided the perfect fluid is a viable model 
for the corresponding multi-particle system consisting of molecules or stars, respectively. 
The equations have been studied accordingly, from mathematical as well as from physical 
viewpoints, see e.g. JH-lEQ. Key problems are existence, uniqueness, (axial) symmetry, 
stability and parametrisation of the solutions, where uniqueness may be understood modulo 
global parameters, like mass and angular momentum. Other useful parametrisations are the 
pressure or the density at the centre or at the the axis, depending on the symmetry, or, as we 
shall see below, the sum of the gravitational potential (p g and the centrifugal potential (p c on 
the surface, called O s . 

In the present work we assume a connected fluid region and focus on investigating the 
problem if the Newton-Euler system admits (only) solutions of finite or infinite extent (i.e., 
compact support of the density function). In fact this question can be regarded as a special case 
of the problem of obtaining a formula for <D S or a bound thereon, for which we give results 
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as well. Normally finite extent of solutions is a prerequisite for the physical significance of a 
model. This always applies to models for single stars. However, for the polytropic EOS 

1 1 n+l 

P = -po" p g " , Po = const. >0 (1) 

n + l 

with index n = 5, all static configurations are infinitely extended. They are known as 
"Plummer's model" for non-rotating globular clusters of stars (see e.g., (H). 

In the static case, systematic treatments of the finiteness problem, which do not rest on the 
assumption of spherical symmetry, can be found in [|5l|6ll3; for stronger results in spherical 
symmetry, cf. (H. Formally, the task consists of determining if the surface value O s of the 
gravitational potential cp g agrees with its value at infinity, which we set to zero. Physically, the 
key ingredient is the continuum form of the virial theorem, and mathematically it is a modified 
Pohozaev-Rellich identity. Moreover, since p g > 0, the maximum principle is available as 
well. In [El|7l the functional 

C p dp' 
Jo Pg(p ) 

was identified as the crucial quantity in the sense that F < for all p guarantees finiteness of 
the solutions unless F = 0, which characterises the polytropes of index 5. Moreover, F > 
for all p but F £ implies that there are no solutions with finite mass. If F changes sign, the 
analysis of the spherically symmetric Newton-Euler system is much more involved and uses 
dynamical systems techniques (see e.g., [9]). 

In the stationary case, results on finiteness of which we are aware are basically 
amendments to theorems on existence (cf. the classic paper OH, the recent account |fTTfl 
and the references therein). All these results require axial symmetry and conditions on the 
EOS which are in many respects more restrictive than ours. In contrast, the present paper 
ignores the problem of existence, which means that we are interested in "a priori" estimates. 
On the other hand, some of our results do not require axial symmetry, which does not hold 
in general for stationary perfect fluids (cf. Sect. 3 below) and we use appropriately adapted 
conditions on the EOS and the LOM. Formally, we compare now the value of the "effective 
potential" (p = 0g+^c at infinity (where we can again set it to zero unless it diverges) and at the 
surface (0 5 ), and the main technical tool is still the virial theorem. The quantity characterising 
rotation which arises in our analysis is 

d 1 

D = x>—d> c + -<p c , (3) 
ox 1 2 

where x' denote Cartesian coordinates. Finiteness of the solution is guaranteed provided D 
and F have the same sign everywhere but do not both vanish identically. Under the additional 
requirement that A(f> c > the case with D and F being non-negative can be excluded by the 
maximum principle, again unless both F and D vanish everywhere. This latter case leads to a 
particular LOM for the n = 5 polytrope, which will be analysed separately. 

For polytopes with index n which rotate according to a "power law" U oc r~ m d/d(p, where 
r is the distance from an axis and m e R, our conditions F < and D < enforcing finiteness 
read n < 5 and m > 5/4, respectively. In particular, a polytropic fluid with n < 5 whose 
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layers rotate at or near their Kepler orbits (m 3/2) must be finite. This is somewhat counter- 
intuitive for the following reason: For such a fluid the gravitational attraction between its 
"layers" will be balanced locally only by the centrifugal force and not by pressure, whence 
the fluid should behave like "dust". However, there is no obvious reason why Keplerian orbits 
of dust (-particles) should not extend to infinity. Of course this "paradox" must disappear 
upon properly taking into account the gravitational interaction between the dust particles or 
fluid layers. However, this reminds of heuristic arguments trying to estimate the amount of 
dark matter necessary to stabilise galaxies and to explain the observed rotation law m « 1. 
Recent analyses of the dark matter problem by different methods (see e.g., |[T2l and IfTBl ) 
indicate that for this purpose much less dark matter is needed than previous approximations 
suggested. In resolving this issue Vlasov-Poisson theory should play a key role and also 
establish a connection to the phenomena in fluid mechanics described here (see e.g., IfPflO . 

While we believe that our results are new, they are rather elementary, and the axially 
symmetric ones might well be contained in the vast astronomical, physical and mathematical 
literature on rotating fluids since Euler's and Newton's time. On the other hand, our exposition 
is motivated by and partially adapted to the corresponding relativistic problem. In fact for 
static perfect fluids in general relativity the quantity F{p) in ([2]) can be replaced by 



to obtain analogous conclusions regarding finiteness and infiniteness as in the Newtonian case 
013. However, since the known relativistic virial theorems |fT51[T6l are not suitable for the 
present purpose the methods are different, and generalisations to the stationary case are not 
straightforward. We intend to present relativistic analogues of some of the results given below 
elsewhere. 

2. Assumptions and Basic results 

This section is divided into five subsections. In Sect. 2.1 we give some basic definitions 
and continue in Sect. 2.2 with writing the Newton-Euler system in a form suitable for our 
purposes. Sect. 2.3 contains a discussion of the EOS p g = p g (p) and the relationships p g ((f>) 
and p((f>). (The symbols were defined in the introduction). The space-time dependence 
of the gravitational variables p g , (p g and the rotational variables p c , <f> c will be discussed in 
Sects. 2.4 and 2.5, respectively, and the latter section also contains the key Lemma 2.5.2 on 
(in)finiteness. 

We denote Cartesian coordinates by either x' (z = 1,2, 3) or by x,y and z. R = 
Vx 2 + x 2 + z 2 , & and ip are spherical polar coordinates, and r = -\Jx 2 + y 2 and (f are cylindrical 
polar coordinates. The volume element on R 3 is denoted by dv. 

2.1. Finiteness versus infiniteness 




(4) 



The vacuum and the fluid regions *V and T are by definition 3-dimensional, open sets where 
p g = and p g ^ 0, respectively. This means that all points ve'V and / e f have open 
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3-neighbourhoods in "V and T . We take the T to be connected but possibly with non-trivial 
topology; in particular *V and S = dT may be disconnected. The boundary S is required to 
be a C 1 submanifold with T lying only on one side of dT. In other words, p g only vanishes 
in vacuum and possibly at the boundary, but not on sets of dimension less than three "inside 
the region occupied by the fluid." The reason for this requirement is to avoid trouble with 
integrating Euler's equation (|9]). We also note that the boundary can be characterised by 
vanishing pressure, cf. Sect. 2.3. 

The following definition introduces a shorthand for our main issue. 

Definition 2.1. The fluid region T is called finite if it stays within a compact subset of R 3 , and 
infinite otherwise. 

In particular, fluids which extend to infinity in at least one direction are called infinite. In 
the axially symmetric case dealt with in Sect. 3 we will distinguish between (in)finite extent 
in axial and equatorial directions. 

2.2. The Newton-Euler system 

We denote by U the velocity of the fluid and by I the integral 

I- f 4^- (5) 

Jo Pg(P') 

(the specific enthalpy), whose existence is assumed for finite p; this is satisfied in particular 
for polytropes p g (p) oc p a when a < 1. Further restrictions on the EOS are discussed in 
Sect. 2.3. The stationary Newton-Euler system can be written as follows 

Pg = PgiP), (6) 

Atffg = 4np g , (7) 

v(p g u) = 0, (8) 

-(£/.v)t/ = V^ + ^ = v(^ + j). (9) 

Pg 

In terms of the centrifugal potential <p c defined up to a constant (which will be specified 
in Sect. 2.5) by 

(f/.V) U = V0 C , (10) 
(|9]) yields the "Bernoulli" equation 

<f> g + <f> c +I(p) = ® s , (11) 
where 5 is a constant. We also introduce a "centrifugal charge density" 

Pc = ldiv[(£/.V)£/], (12) 
in terms of which Eq. ([9]) yields 

A0 C = 4np c . (13) 
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Remarks. 

(i) The existence theorems in the axially symmetric case show that the motion of the fluid 
has to be specified somehow in the Newton-Euler system. A simple way of doing 
so is to prescribe the velocity in terms of position U = U(x') or to prescribe the 
centrufugal potential cp c (x') which we do below (without restriction to axial symmetry). 
An alternative is to specify the angular momentum per unit mass (see e.g., IfTOlO . 

(ii) The above definitions of U, (p c and p c , and the equations containing them, are understood 
to hold only in f. It will however be convenient to extend these quantities, and some 
relations between them, to *V. In Sect. 2.5 we will discuss four alternative extensions 
which will be used in the following sections. 

(iii) In work focusing on existence, the aim is to specify only the EOS and the velocity or the 
specific angular momentum as functions in space, and to get information on the spatial 
behaviour of all variables. In this respect the present work has the same scope, although 
existence is not the issue here. However, except for crucial conditions on the EOS and the 
LOM, we will in the sequel also have to make differentiability and falloff requirements 
for the space-time dependence of our functions. 

2.3. The Equation of state and the effective potential 

The "effective potential" is defined by cp = <p g + <p c ; its gradient is called "effective gravity" in 
Sect. 3.2.1 of ED. 

Our first Lemma serves mainly to list the assumptions on the equation of state required 
later. The proof is an easy consequence of the Bernoulli equation ( fTT| ). 

Lemma 2.3.1. We assume that p g (p) is piecewise continuous, < p g (p) < oo, and that the 
integral I exists for finite p. Then in T the effective potential (p(p) is C°, piecewise C 1 , and 
strictly monotonic; the same applies to the inverse p = p((p - 5 ) = p((f>), and the density p g 
is also a C° and piecewise C 1 function of (p that satisfies p g = dp((p)/d(p. 

Remark. The Lemma implies that the surface p = is an equipotential surface of (p, and 
4>{p = 0) takes the value 5 . If we allowed for disconnected fluid regions, a consistent 
definition of <f> on R 3 (cf. Sect. 2.5) would imply different constants <D 5 on each component 
in general. 

The next Lemma (which is known, see e.g., IfTTll ) contains a stronger assumption on the 
EOS in a neighbourhood of p = 0, which will also be made in Proposition 4.2. 

Lemma 2.3.2. In addition to the requirements of Lemma 2.3.1, we assume that p g (p) is C° in 
[0, 6) and C 1 in (0, 8) for some 6 > 0. Then lim p ^ plp g ip) = 0. 

Proof. The result is obvious if p g (0) ± 0. If p g (0) = we first note that near p = the inverse 
p = pipg) exists, and dp/dp g > 0. We can thus replace the assertion by lim p ^ pip)lp g = 
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(dropping the subscript g on p g when the latter is a sub- or superscript itself). We obtain 



oo > lim 



. r p dp' r Psip) 

im I = lim I 

^°Je Pg(P') ^Jp(e) 



Pg(p') 

dp(p g ) 



P 

> lim 1 - 6 ' lim 

p->o Jp § ^°J P (6) P' g p->° dp g 



dp'dp' g 

— > hm 

^Jpie) dp' g p' g 

*v>df/ K dp{p g ) 

= hm . oo 



inf 



dpiPg) 



\Pg{e),p g {p)\ dpg 



r p(p) d A\ > 

J pie) P'g J 



(14) 



Hence = lim p ^ dp(p g )/dp g = lim p ^ p(Pg)/p g by de l'Hospital's rule. 



2.4. The gravitational variables 

We work in the weighted Sobolev spaces W g ' p (\ < p e R, 6 e R, k e No), based on the 
weighted Lebesgue norms 




for measurable functions u 6 L^ c (R 3 j, where cr = (1 + 7? 2 ) 1/2 . This is Bartnrk's index 
convention, cf. lfT9l . see also remark (ii) below. 

We always require that the fluid has finite mass m = p g dv. The following Lemma is 
standard; we use [fT8l[T9ll for the inversion of the Laplacian in ([7]), and the strong maximum 
principle, Thm. 9.6 of flU. 

2 2 2 J~~ L 

Lemma 2.4. Let p g 6 W^_3_q. 5 < a < 1. Then there is a unique solution <p g 6 W /oc of 071) with 
iffg = (p g - m/ a 6 Wj^. Moreover, g < 0. 



Remarks. 

(i) Weighted Sobolev spaces have the important property that a rather slow falloff of the 
density function in one or two directions is admitted as long as it is compensated by 
sufficiently fast falloff in the other direction(s). In particular, an axially symmetric 
disk of finite thickness with a measurable density function p g 6 W ' 2 = L 2 and falloff 
Pg = 0(l/r 2+e ), 6 > 0, satisfies the requirement of Lemma 2.4 and of the subsequent 
results. 

(ii) With the conventions of Bartnik ffT9l used above, the index yS for / e W^' p is related 
to the growth of / at infinity; in particular / = o(R^) provided kp > 3 (throughout 
the paper, falloff conditions are always understood for large r or R). In fact for the 
function \f/ g introduced in Lemma 2.4 it follows that tf/ g = o(i?~ l e ). However, getting the 
corresponding falloff for the derivatives, namely V,i/r g = o(R~ 2 ~ £ ) would require p > 3. 
While such first derivatives do occur in Theorem 4.2 below, less precise information on 
their falloff, which follows from p = 2, will suffice. 



2.5. The rotational variables 

As already mentioned in the introduction, the strategy of our finiteness argument is to compare 
the value of the effective potential = (p g + <p c at the surface with its value at infinity. On the 
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other hand, integrating Euler's equation (|9]) defines <p c only in T, and in fact only up to a 
constant. 

We therefore extend now (p c to ( V. Below we consider four alternative definitions of 
the rotational variables valid on R 3 , labelled A, B<, B> and C. In extensions B<, B> and 
C we follow standard practice and prescribe <p c (x l ) a priori on R 3 , irrespective of the fluid 
distribution. Of course we can also prescribe U(x l ) instead of (p c (x l ), together with a constant 
in the resulting <p c . On the other hand, extension A is motivated by the analogy between the 
Poisson equations ([7]) and ( 13 ), which suggests treating p c as a "source" for <p c . 



While A, B< and B > are compatible with any symmetry, definition C refers only to 
a cylindrically symmetric velocity distribution (f> c (r). All definitions are in fact not only 
extensions from T to <V, but involve also extra conditions on the falloff of the velocity at 
infinity if the fluid spreads there (cf. the remarks after Definition 2.5 for details). 

Definition 2.5. For a solution of ([6 -11) with I finite, p g and (f> g as in Lemma 2.4, and 



1 < q < oo, < e < 1, we assume one of A, B or C: 



A: p c 6 W_2_ in T , p c = in "V and <p c — » at infinity. 

B: <p c {x l ) e W 2,q given on R 3 such that, for all radial unit vectors n there is a unique limit 

lim <f> c (R ft) = = <f>U&, <P) (16) 

R — ><x> 

and either 

B<: sup^0oo(^,^) = 0, or 
B>: mf^fiW&TP) = 0. 
C: cf> c (r) eC'(R\ {0}) given (and possibly divergent at 0), and lirn,--^ <p c {r) = 0. 

Remarks. 

(i) While condition A is somewhat alien to the Newtonian case, it is mainly motivated by 
Relativity. There the norm and twist potentials of the stationary Killing field satisfy 
elliptic equations, while on the other hand counterparts of the coordinate conditions B 
and C will hardly make good sense. 



(ii) We note that cylindrically symmetric potentials <p c (r) do not satisfy (16) unless (p c {r) = 
const., since the limit lim^oo 4> c (r) = lim^o, (f> c (r) = c (r) in the axis direction depends 
manifestly on r, while there should be a unique 0co(# = 0,0). On the other hand, 
condition B is meaningful in particular for "almost spherical" velocity distributions. 

(iii) Since <p c was defined in ( [TT] ) only up to a constant, condition B amounts to requiring 
( 16), and that either sup#^ 4>oo{$, <f) < oo, or inf^ 0oo(#, (p) > -oo, while C includes the 



requirement that lim r ^oo <p c (r) > -oo. This latter requirement can in fact be removed, so 
that (f> c (r) is allowed to diverge both at the axis and at infinity. However, in order not to 
overload the subsequent definitions and results, this option will be considered only in the 
final remark of Sect. 3.2 and in examples. 

(iv) The falloff conditions B< and B> and C are in some sense less restrictive than A. However 
in Theorem 4.2 the former requirements need to be explicitly supplemented by a falloff 
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condition on the pressure which, on the other hand, follows from A automatically by 
virtue of the Bernoulli equation. 

(v) None of the definitions A, B<, B> or C extends the validity of the Bernoulli equation 
from T to the vacuum region *V in general. Such an extension could simply be afforded 
by setting <p c = -<p g + 5 in "V, but it seems to be of little use. Note in particular that 
such a 4> c would not necessarily be C 1 at the surface. 

(vi) In principle, we could also consider axially symmetric velocity distributions of the 
form (f> c (r,z). However, the Poincare-Wavre theorem, Lemma 3.1, then already implies 
cylindrical symmetry in T , so an extension to "V with the same symmetry is the natural 
choice. 

The following Lemma is analogous to Lemma 2.4 but, compared to the gravitational 
variables, the falloffs are slower here, and in contrast to the natural condition p g > 0, the 
assumption p c > is highly restrictive. 

Lemma 2.5.1. If p c e W_f_ , or > 0, then <f> c can be chosen such that (p c e W^. Moreover, if 
p c > 0, then cf) c < 0, and therefore also 5 < 0. 

We now have the following easy to prove, but important Lemma. 

Lemma 2.5.2. Under the requirements and with the labelling of Definition 2.5 the following 
holds: 

A: If ® s # then the fluid is finite. Moreover, if p c > the fluid is finite if <J> 5 < and 
infinite if 5 = 0. 

B: If 0co(#, ip) does not agree with <D S for some ip), the fluid is finite "in the direction 
ip)" or more precisely, it does not intersect any 2-sphere of sufficiently large radius 
R at the points {R,-&,ip) e S^. In particular, for extensions B<, B> the fluid is finite if 
O s > 0, 5 <0 respectively. 

C: If <D S < 0, the fluid is finite (in all directions), and if 5 > the fluid is finite in every 
direction except possibly in the axial one. Moreover, if <D S = 0, and the fluid extends to 
infinity in the axis direction at some radius r { , it is static for all r>r x . 

Proof. The proofs of parts A and B follow easily from the definitions, and the second part of 
case A from the maximum principle. The proof of C is postponed to Section 3.2. 

Remark. We recall from a previous remark that a cylindrically symmetric cf) c (r) considered 



under C in Lemma 2.5.2 is not compatible with B, in particular not with ( 16 ). This accounts 



for differences in the conclusions and the proofs in the corresponding parts of Lemma 2.5.2. 



3. Results with Symmetry 

Assumptions and results on symmetry of barotropic perfect fluids require a careful discussion 
(cf. lETI where the relativistic case is included as well). 

It is important to distinguish between the symmetry of the velocity field U(x l ) and its 
trajectories on the one hand, and the symmetry of the whole configuration on the other 
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hand. In fact, a well known example important in theory are non-axially symmetric fluid 
trajectories winding around an axially symmetric torus (cf. Q]| in the Newtonian case and Il22l 
in Relativity). Fluids with non-axially symmetric configuration have been discussed; cf. |j2T| 
and Sect. 2.8.3 of Q. On the other hand, known existence proofs of reasonable generality all 
assume axially symmetric velocities (cf. IPTO and the references therein). 

In the Newton-Euler system ([6]-[9]) the gravitational and velocity variables can enjoy 
different symmetries, as long as the coupling between the potentials via the Bernoulli equation 



(11) is respected. In this work we do not assume any symmetry of the matter variables p g 
and (p g . As to the rotational variables we first state in Sect. 3.1 the definitions of axial and 
cylindrical symmetry. Then we continue with a known result (the "Poincare-Wavre theorem") 
relating them, and give a simple application. In Sect. 3.2 we examine systematically the 
results arising from the integrals of the cylindrically symmetric Bernoulli equation via the 
limits r — > oo and z — » °°. In particular, Sect. 3.2 also contains the proof of Lemma 2.5.2.C. 

3.1. The Poincare-Wavre theorem 

Definition 3.1. This refers to solutions to the Newton-Euler system ([6]-[9]): 

(i) A solution has axially symmetric velocity if the velocity tJ is proportional to the axial 
Killing vector, and if \U\ (and hence also <p c ) are rotation invariant. 

(ii) An axially symmetric solution has cylindrically symmetric velocity if U, (and hence also 
(p c ) are invariant under translations along an axis. 

Lemma 3.1. For a solution of Newton-Euler system (|6]-[9j) with axially symmetric centrifugal 
potential (p c (r,z), the latter is in fact cylindrically symmetric, viz. 

cf> c (r,z) = (f>c(r), = u*r)—, (17) 

dip 

and it holds 

(f> c (r) = - f co 2 (r')r'dr', r > r { = const. (18) 

Proof. The proof is obvious from Euler's equation which reads 
d(f> c (r,z) d<p c (r,z) 2 

— 7 = 0, — = -co(r,z)r (19) 

oz or 

under the stated assumptions. 

Here is a simple application of the previous Lemma. 

Proposition 3.1. We assume a cylindrically symmetric velocity distribution and an EOS as in 
Lemma 2.3.1. If the fluid extends to infinity in the axis direction at two different radii r\ > 
and r 2 > r\, the intermediate region J{ = {(r,z) with r\ < r < r 2 and - oo < z < oo} must be 
static, i.e., cp c = const, in Si. 

Proof. From the Bernoulli equation we obtain 

<f>M) = lim (p c (n) = lim (<Dj - ^ - j) = O s , (20) 

7— > CO 7— *CO * ' 



(In)finite extent of stationary perfect fluids in Newtonian theory 



10 



and analogously (p c (r 2 ) = 5 , which contradicts ( 18 ), unless to = 0, and the solution is static. 

Remark. This result implies in particular that solutions with axially symmetric velocity 
distribution which fill the whole space must be static and hence spherically symmetric. 

3.2. Finiteness and <D 5 

In this section we discuss the finiteness question under the assumption that we know the sign 
or the vanishing of O5; in particular we will prove Lemma 2.5.2.C. On the other hand, in 
Sect. 4 the virial theorem will be employed to obtain the required information on <D S from the 
EOS and the LOM, independently of symmetry assumptions. 

Finiteness in the radial direction, and the behaviour near the axis, can be obtained rather 
easily from a qualitative discussion of the form of the potential. The key features of the 
potentials to be kept in mind here are: 

(i) <f> g < 0. 



(ii) 4> c {r) behaves as in Definition 2.5.C and is monotonically decreasing with r from ( 18 ). 

(iii) Inside the fluid region cp = tp g + (p c <<& s , and tp /" ^s* as one approaches a locus of zero 
pressure (irrespective of its location in space). 

Figures l.(a-d) show qualitatively possible forms of the potentials together with the 
respective fluid fillings (shaded), differing in their range with respect to infinity and the axis. 
For positive potentials tp as in Fig. l.(a) no solution exists, irrespective of the value of <D 5 , 
while the other figures show ranges for tp which are a priori admissible but still without 
guarantee of existence of a solution. If the fluid is not axially symmetric, the curves are 
understood to be sections of higher dimensional level surfaces of the potentials. 

While the behaviour of the fluid in the axis direction as included in Lemma 2.5.2.C is 
not seen from the above diagrams, it is a simple consequence of the Bernoulli equation. The 
Lemma is now proven as follows. 

Proof of Lemma 2.5.2.C. By contradiction: Assume that the fluid is infinite in the r direction. 



Then (11) and Definition 2.5.C yield 

5 = lim (cp + I) = 0, (21) 



which contradicts the behaviour of O s , as claimed for finiteness in the r direction. Here ( |2T| ) 
in fact holds for the limit in any direction not parallel to the axis. The final assertion of the 



Lemma is obtained from (11), ( 18 ) and Definition 2.5.C. as follows: Assuming that the fluid 



extends to infinity in the z direction at r we get 

/~»oo 

5 = lim (<p + 1) = <p c (r) = I to 2 (r')r'dr'. (22) 
and hence O s = implies to = 0. 

Remark. As remarked after Definition 2.5, the condition in part C that lim^oo <p c {r) > -00 
can be removed. In fact, if lim,.^ (p c (r) = -00, the statement of Lemma 2.5.2.C remains 
correct and the proof is trivially adapted. 




Figure 1. Different a priori possible forms of the potentials <p g , <p c and <p. 



4. General Results 

In Sect. 4.1 we give a result which, under rather restrictive conditions, forbids shells or 
"hollow" bodies. We continue in Sect. 4.2 with our main finiteness theorem, which rests 
on the suitably adapted virial theorem, and we also obtain a more general bound on 5 . 



4.1. A "no-shell" result 



Proposition 4.1. We consider a solution of ([6- 13 ) and assume part A of Definition 2.5, with 
p c > 0. Then there cannot be any vacuum region not connected to infinity, i.e., the fluid T is 
not a "shell", and it is not "hollow". 
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Proof. Assuming the contrary, T has an inner and an outer boundary, the latter either 
separating T from infinity or located at infinity, where <p takes on the value O s . We first 



apply the strong maximum principle, Thm. 9.6 of 11201 . to ( 13 ) on the region % consisting of 
T and the "inner vacuum," but not the vacuum component connected to infinity. This implies 
that (p = O s is the maximum taken at the boundary of % (the outer boundary of T). But this 
maximum is also taken on at interior points of K , namely the inner boundary of f, and hence 
the maximum principle implies that </> is constant on K. The latter, however, is easily excluded 
from the assumptions. 

Remark. Clearly, this result does not exclude toroidal rings of finite thickness, with or without 
central body. 



4.2. The finiteness theorem 

Our main result is now obtained by combining various Lemmas of the preceding sections. 



The quantities F and D were defined in the Introduction; the main conditions (23) and (24) 
are discussed in remarks after the following Theorem and in Sect. 5. 



Theorem 4.2. Assume we are given a solution of ([6 -11) such that the EOS ([6]) satisfies 
the requirements of Lemma 2.3.1, and that p g e W\_ , e > as in Lemma 2.4. Moreover, we 
require that (p c or p c should satisfy one of A, B<, B> or C in Definition 2.5 and that there holds 
in case 



A: one of 



F(p) > Vp and D(x) > Vx, 
F(p)<0 Vp and D(x) < Vx, 



(23) 
(24) 



B<: p e W x * and (1231), 



4-6 

1,1 



B>: pe W^i and (24), 



-4-e 

C: pe W l _'l £ and {24]). 



Then either O s ^ and the fluid is finite, or F = 0, D = and 5 = 0. 



Corollary 4.2.1. In case C above, assume that (23) holds instead of (24). Then the fluid is 



finite in all directions except possibly in the axial one. 



Corollary 4.2.2. In case A above, assume that there hold (24) and p c > 0. Then either <D S < 
and the fluid is finite, or F = 0, D = and <D 5 = 0, and the fluid is infinite. 

Proof. We use the following modified version of the Pohozaev-Rellich identity []24ll25l : Let 
|=z'|ibea dilation, i.e., V (i ^ ;) = ^(V,^- + V y £) = where gij denotes the components of 
the 3-metric, and V, is the covariant derivative with respect to Then a simple calculation 
(compare [|5ll3) shows that 



= 2n[p g (<f> g -2?V i( f> t ) + 6p}. (25) 
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We first note that p e W[^_ e , which in case A follows from Euler's equation (9). Next, by 
Gauss's law, the left hand side can be written as a surface term over a ball S of radius R. We 
now insert the forms of <p g , V0 g , obtained in Lemma 2.4, and use the fact that ifr g and V,-^ g 
defined in this Lemma "fall off faster" than the leading terms, to show that the surface integral 
vanishes as R — > oo. This follows straightforwardly via Cauchy-Schwarz estimates and by 
applying the "trace theorem" (see e.g., Sect. 5.5 of Il23~l0 : 

ll/laslb (as) <C||/||^ (<S) (26) 

(which holds for any function / e W 1,q (S), q > 1, and its extension f\ ds to dS), to the 
functions tf/ g , Vflffg and p. In the limit R — » oo we are left with 



= In J 3 [p g (<p g - 2^V,-0 C ) + 6p] dv. 



Inserting (11) into the above expression gives 

mO s = J [F(x) + 2p g D(x)] dv. 
The theorem is now obvious from the requirements and from Lemma 2.5.2. 



(27) 



(28) 



Remarks. 

(i) Note that the conclusion does not require falloff conditions for F(x), and only the mild 
falloff conditions on 4> c {x j ) or p c from Definition 2.5. Even without such conditions, Eq. 



(28 ) and finiteness of mO s imply that the integral on the right exists. 

(ii) In the static case existence of finite, spherically symmetric solutions is known for 
poly tropes with index n < 5, which corresponds to F < 0. On the other hand, if F > 
solutions with finite mass do not exist except for the poly tropes of index n = 5. This 



suggests that the "realistic" range for rotating fluids is given by conditions (24). In fact 
the requirements of existence results for axisymmetric rotating fluids of which we are 



aware ( IfTTTl and the references therein) fall in this range, whereas (23) might not allow 
for solutions at all. 



(iii) If p c is at least weakly differentiable, the second term on the right in ( 28 ) can be rewritten 



as 



I p g D(x)dv= I <j> g 
Jr 3 Jr 3 



5 



dv + surface terms, 



(29) 



and, since (p g < 0, the sign of the expression in the bracket determines (in)finite extent 
of the fluid region in the same manner as D(x'). The surface terms vanish in case A but 
have to be handled with care in cases B<, B> and C. 

(iv) To show the relation with the classical virial theorem (cf. Sect. 1.4 of [fT5l or Sect. 2.8.1 
of 0) we use ([25]) and definition (TTOj) to obtain ([30]) below. To get pT]) which 



is the sum of the potential energy E pot = 1/2 J_ 3 p g <p g dv, the (bulk) kinetic energy 
Eun = 1/2 L 3 p g \U\ 2 dv, and the thermal energy (kinetic energy of the thermal motion) 
E therm = 3/2 pdv, one has to remove the second term in (30) by partial integration 
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using the continuity equation (|8J), and assume that the velocity U falls off suitably in 
order for the surface terms to vanish. 



= 2n f [p g [<f> g - 2U j Vj (eUt) + 2\U\ 2 ] + 6p) dv 



(30) 
(31) 



(v) In the static case 5 is the gravitational potential at the surface, and it is related to the 
observed redshift. For a rotating, extended object, the redshift on the rim clearly depends 
on the velocity of the rotation as well, but does not seem to bear any obvious relation to 

o 5 . 



(vi) Equation ( 28 ) is a formula for <D S , but as such it requires a knowledge of F(x), as opposed 
to the mere sign condition on F(p) used before. An estimate for <D S in terms of F(p) and 
D(x) is the following. 



Proposition 4.2. Under the requirements of Lemma 2.3.2, 

F(p) 



|®sl < sup 

r 



Pg 



+ 2 sup \D(x) 

r 



(32) 



Proof. From (28), 

m|0 5 | 



<-L 



< m sup 

r 



Pg sup 
F(p) 



F(x) 



Pg 



+ 2D(x) 



dv < m sup 

r 



F(Pg) 



Pg 



+ 2D(x) 



Pg 



+ 2m sup \D(x)\ 

r 



(33) 



and the supremum of F(p)/p g (p) exists, since from Lemma 2.3.2 lim^o F(p)/p g (p) = 
lim p ^ plPgip) = . 

5. Examples 

We examine here the limiting case of Theorem 4.2, and conclude with discussing "power law 
rotations." Explicit examples are only available in the cylindrically symmetric case. 



5.1. The limiting case 

The limiting case of Theorem 4.2 is F(p) = 0, D = and <D S = 0. The first condition yields 
the 1 -parameter family of poly tropic EOS ([T]) with index n = 5. In the static case, the resulting 
PDE 



A(p g = -4np (f> 5 g 

has for each p a the well-known (cf. e.g. flU) family of solutions 

m 

= j 

^JL pofn A +R 2 

parametrised by the mass m. All these solutions extend to infinity. 



(34) 
(35) 
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In the stationary case, if p c is differentiable, it follows from = D = x l d^ c + \(p c that 
x'dip c + |p c = 0. This yields that </> c and p c are homogeneous functions of degree -1/2 and 
-5/2, respectively. That is to say, these functions have the form 

0c = z~*<r(~,^), Pc = z^r{^,^y z±0 (36) 

for some arbitrary (but mutually related) functions cr(x,y) and r(x,y). To determine the 
solution, we have to solve 

Acf> = A(cf> g + (p c ) = -4n Po cf> 5 +4jtz-*t^J} (37) 

in T , and A0 = in <V (if present). We remark that (37 ) is scale invariant under 

0(jc') — > J -(/> | - ) , V A = const. > 0. (38) 



Thus the motion determined by D = could be called "scale invariantly rotating 
poly trope of index 5." 

Scale invariance is sometimes useful for getting information about (non-)existence of 
solutions, in particular in combination with the scaling behaviour of the energy functional. 
We do not go into details here. 

Under restriction to axially symmetric rotation laws (about the the z axis) (36 ) becomes 

c = r*arQ, p c =z~^fi^j, z±0 (39) 

for some arbitrary (again related) functions a(r) and yS(r). However, by Lemma 3.1, the 
rotation law must be cylindrically symmetric. Choosing a and f3 in (39) appropriately, we 
have 

C 2 C 2 

0,(r) = 2—, p c (r) = T , C = const, in T. (40) 

yr 87rn 

We obtain the following behaviour for the "scale invariant rotation" o» s ,(r), whose falloff 
interestingly lies between the Kepler angular velocity co Q (actually Copernicus ll26ll . as this 
is for circular orbits) and the observed galaxy rotation curves u)@: 

c c c 

co Q = — , co si = — , a>@ . (41) 

Yi f4 r 

Regarding finiteness, we note that co and 4> c are singular on the axis, whence Theorem 
4.2.A is not applicable; in particular the maximum of </> is obviously taken at the axis. 
However, d40j) implies that (f> = (p c + (f> g = 2C 2 / - m/R + o(l/R l+£ ) > for sufficiently 
large r. Hence if a solution exists, the form of must be as in Fig. l.(c), while the form of 
Fig. l.(b) is excluded at least in the cylindrically symmetric case considered here. Thus the 
fluid is either finite (with a hole near the axis) or it is infinite in the z-direction (respecting 
Proposition 3.1). The former case would probably give a toroidal configuration. 

We can now discuss the consequences of the above arguments for "power law rotations." 
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5.2. Power law rotation 



Proposition 5.2. Assume ([6-11 ) have a solution which satisfies the conditions of Theorem 
4.2.C with the rotation law 

* c(r)= 2(k- C V)r« k -» ^ W(r) = ^' KCeR - (42) 
Then depending on the value of k the fluid has the following properties regarding infinity and 
the axis: 

(i) k < |: It is finite in the radial direction. 

(ii) k > 1: It has a "hole" near the axis. 

(hi) k > |, F < and p e W_4_ e : It * s either finite (with a hole near the axis), or infinite in 
the axis direction (respecting Proposition 3.1). 

Moreover, if the solution is finite in the axis direction in case (iii), the conclusion (finiteness) 
holds without the falloff condition on the pressure. 

Proof. For k < 1 statement (i) follows form the fact that (p c diverges at infinity. For 1 < k < 3/2 
the same conclusion can be inferred from the argument used in Sec. 5.1 to show that for the 
"scale invariantly" rotating polytrope (p = 4> c + <p g > for sufficiently large r. Conclusion (ii) 
is obvious from the fact that C diverges at the axis. To see (iii) we calculate from Definition 

and use Theorem 4.2.C for the generic case, and the discussion of Sect. 5.1 for the limiting 
case with <D 5 = 0. For the final statement, we note that ([9]) with Lemma 2.4 and k > 5/4 
imply that p = 0(r~ 7/2 ). This is sufficient in order for the surface term in the integral of (25 ) 
to vanish as R — > oo, provided the fluid is finite in the axis direction. 
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